Continuous percolation phase transitions of random networks under a generalized 

Achlioptas process 
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Using the finite-size scaling, we have investigated the percolation phase transitions of evolving 
random networks under a generalized Achlioptas process (GAP). During this GAP, the edge with 
minimum product of two connecting cluster sizes is taken with a probability p from two randomly 
chosen edges. This model becomes the Erdos-Renyi network at p = 0.5 and the random network 
under the Achlioptas process at p = 1. Using both the fixed point of S2/S1 and the straight line 
of lnsi, where Si and S2 are the reduced sizes of the largest and the second largest cluster, we 
demonstrate that the phase transitions of this model are continuous for 0.5 < p < 1. From the 
slopes of lnsi and ln(s2/si)' at the critical point we get the critical exponents ft and v, which 
depend on p. Therefore the universality class of this model should be characterized by p also. 
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two connecting cluster sizes is taken with a probability 
p |22j. At p = 0.5, this model becomes the Erdos-Renyi 
(ER) network [Hj]. At p = 1, the GAP is equivalent 
to the Achlioptas process which suppress the appearance 
of larger clusters. With the increase of p from 0.5 to 1, 
this suppression is switched on gradually. After inves- 
tigating the percolation phase transitions at different p, 
we can get more understanding about these phase tran- 
sitions and make a more reliable conclusion about the 
character of the percolation phase transition at p = 1. 

In our Monte carlo simulations, we begin with N iso- 
lated nodes and then connect them with the edges added 
through the GAP. The network obtained is character- 
ized by N, the number of edges N r and the probabil- 
ity parameter p of GAP. For p < 1, we take N = 
(4, 8, 16, 32, 64, 128) x 10 4 and run 500,000 steps in 
each simulation. At p = 1, we choose much larger system 
sizes N = (512, 1024, 2048) x 10 4 to reach the asymp- 
totic region of finite-size effects and run 100, 000 steps in 
each simulation. 

For the cluster ranked R and with size S R (N r , N; p), 
we define its reduced size as 



The modern theory of complex networks [l|-|3( has 
opened new perspectives in the study of complex sys- 
tems in nature, social and economic systems, technical 
infrastructures and many other fields. The macroscopic 
properties of complex networks emerge from the inter- 
actions among individual constituents. The percolation 
phase transition of complex networks is an interesting 
macroscopic property and can be studied with the per- 
colation theory in statistical physics It has been 
pointed out that the random network undergoes a con- 
tinuous percolation phase transition during the random 
process |5j . The critical phenomena of complex networks 
are reviewed in Rcf . [6] . 

However, Achlioptas et al.Q reported that the perco- 
lation phase transition of random networks becomes dis- 
continuous under the Achlioptas process (AP). During 
this AP, two unoccupied edges are chosen randomly and 
the edge with minimum product of the connecting cluster 
sizes is taken as the next occupied bond. This discontin- 
uous phase transition in networks is called the explosive 
percolation Q. The AP was applied later to other net- 
works and some similar rules were also introduced 
[§-[l5j]. Using the finite-size scaling, the explosive per- 
colation phase transition was also investigated [ill [l3] ■ 
These works support that the explosive percolation is 
discontinuous. Recently, Costa et al.fljl] showed that the 
explosive percolation transition is actually a continuous 
phase transition with a uniquely small critical exponent. 
Riordan and Warnke [2(| prove mathematically that all 
Achlioptas processes have continuous phase transitions. 

In this letter we investigate the percolation phase tran- 
sitions of random networks under a generalized Achliop- 
tas process (GAP). In the GAP, two unoccupied edges are 
chosen randomly and the edge with minimum product of 



s R (r, N; p) = S R (N r , N; p)/N 



(1) 



where r = N r /N is the reduced edge number. The 
reduced size of the largest cluster si(r,N;p) at N — 
1.28 x 10 6 is shown in Fig.[TJ si becomes finite for r > r c 
and there is a percolation phase transition. The tran- 
sition point r c increases with p and the corresponding 
phase transition becomes sharper. 

If the percolation phase transition is continuous, the 
reduced sizes sr(v, N;p) should follow a finite-size scal- 
ing form [illlil 



s R (r,N;p) = N-^J R (tN 1 ^;p), 



(2) 
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Reduced size of the largest cluster at N — 1.28 x 



where t = (r — r c )/r c and v is the critical exponent of the 
correlation length £ = £o '• The scaling form in Eq. 
([2]) is valid in the asymptotic critical region with \t\ <C 1 
and N 3> 1. Outside this region, additional correction 
terms should be taken into account. 
From Eq. ©, we can obtain 



s 2 /s 1 =s 2 (tN 1 ^;p)/s 1 (tN 1/u ; P ) = U(tN^ v ;p) (3) 



and 



lnsi(r,iV;p) = -(i / v\nN + \ns 1 {tN 1/u -p). (4) 



At a critical point, S2/si\ r =r c = U(0;p) is a fixed point 
versus N and lnsi(r c , 2V;p) = — /3/VlniV + lnsi(0;p) is 
a straight line versus IniV. Using these properties, the 
critical point of complex network can be determined both 
from the fixed point of 82/si and the straight line of 
lnsi versus In N . The critical exponent ratio j3/i> can be 
obtained from the slope of lnsi at r c . 

For p = 0.5, our model is equivalent to the ER model 
25j. Its continuous percolation phase transition can be 
studied with S2/ s% and In s\. In Fig.[2][A), S2/ s% is shown 
as a function of r at different N. The different curves of 
S2/S1 for different N have a fixed point at r\ = 0.5000(4). 
In Fig. HJB), lnsi is shown versus In A for different r. 
The curvature of In si is negative at r = 0.4996 and be- 
comes positive at r — 0.5004. At r — 0.5, lnsi is a 
straight line with slope -0.33(1). We get r 2 c = 0.5000(4) 
in consistence with r\. Our results of r c and /3/v agree 
with the exact results r c = 0.5 and P/v = 1/3 of the ER 
model [H. 

In Fig. [3l we show s^/si and lnsi of random network 
under the GAP at p = 0.8. From the fixed point of 
s 2 /si, we get the critical point r\ = 0.7252 ±0.0001. The 
curvature of lnsi is negative at r — 0.7250 and becomes 
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Figure 2: (A) s 2 /si has a fixed point at rl = 0.5000 ± 0.0004; 
(B) Log-log plot of the reduced size si versus N. Around 
r'i = 0.5000, the curvature of lnsi changes its sign. 



positive at r = 0.7254. So r\ = 0.7252 ± 0.0002, which is 
consistent with r\. 

For p — 1.0, our model becomes the PR model of 
Ref.[7J. s 2 j s\ at different N are shown in Fig. BJA) an d 
there is a fixed-point at = 0.88844(2). The curvature 
of In si is negative at r = 0.88840 and becomes positive at 
r = 0.88850. We obtain r 2 c = 0.88845(5) consistent with 
t\. From the slope of lnsi at r 2 , we get ft/v = 0.04(1). 
Our r c and f}/v at p = 1 agree well with Ref.[26]. 

To determine the critical exponent v of correlation 
length, we study the derivative 



dr 



= N 1/v r- 1 U'{tN 1/v -p) 



At a critical point, we have 

ln(s 2 /si)'|r=r c = [1/v) lnN + lnU'(0;p)-lnr c 



(5) 



(6) 



which is a straight line with slope 1/v. From our Monte 
Carlo data, we can calculate 1/v at different p, which are 
summarized in Table I. 
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Figure 3: (A)s 2 /s 1 has a fixed point at r\ = 0.7252 ± 0.0001; 
(B)Log-log plot ol the reduced size si versus N at different r. 
Around r 2 = 0.7252, the curvature of In si changes its sign. 



Figure 4: (A) s 2 /si has a fixed point at r\ = 0.88844 ± 
0.00002. (B) Log-log plot of the reduced size si versus N. At 
r'i = 0.88845, the slope of lnsi versus In AT gives fi/v « 0.04. 



Table I: Critical reduced edge number r c , critical exponent 
ratio fifv and inverse of critical exponent 1/u. We obtain r\ 
from S2/S1 and r\ from lnsi- 



p 


r 1 

1 C 


r 2 


p/u 


1/v 


0.5 


0.5000(4) 


0.5000(4) 


0.33(1) 


0.33(1) 


0.6 


0.5599(1) 


0.5599(3) 


0.32(1) 


0.33(1) 


0.7 


0.6349(2) 


0.6349(3) 


0.32(1) 


0.34(1) 


0.8 


0.7252(1) 


0.7252(2) 


0.33(1) 


0.35(1) 


0.9 


0.8216(1) 


0.8216(3) 


0.31(1) 


0.38(1) 


0.95 


0.8603(2) 


0.8604(2) 


0.15(1) 


0.46(1) 


1.0 


0.88844(2) 


0.88845(5) 


0.04(1) 


0.50(1) 



Using r c and \ jv given in Table I , we can define a 
finite-size scaling variable tN 1 '" at each p. With this 
scaling variable, our Monte Carlo simulation data of 
si I Si(r, N;p) at a given p and different N collapse into a 
finite-size scaling function {/(tN 1 ^-^), which are shown 
in Fig. [S] for p = 0.5 and p = 1. The p-dependence of 
U(0;p) is demonstrated in Fig. EI 



In summary, we have studied the percolation phase 
transitions of evolving random networks under a gener- 
alized Achlioptas process that the edge with minimum 
product of two connecting cluster sizes is taken with a 
probability p from two randomly chosen edges. This 
model becomes the ER network [25| at p — 0.5 and 
the random network under the Achlioptas process [7| at 
p = 1. Using the finite-size scaling, the percolation phase 
transitions of this model are studied. From the finite-size 
scaling forms of the largest cluster size s\ (r, N; p) and the 
second largest cluster size S2(r,N;p), the fixed point of 
S2/S1 and the straight line of lnsi versus lnA^ can be 
used to determine the critical points. It has been found 
that the critical points r c determined from si/si and 
In si are consistent and increase with p. From the slopes 
of lnsi and ln(s2/si)' at r c , we can obtain the critical 
exponent ratio j3/v and 1/v respectively, fijv decreases 
from 0.33(1) at p = 0.5 to 0.04(1) at p = 1. 1/v increases 
from 0.33(1) at p = 0.5 to 0.50(1) p= 1. With the scaling 
variable tN x / v , the Monte Carlo data of S2I s\(r, N;p) at 
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Figure 5: Finite-size scaling function U(tN 1,lv ;p) at (A) p 
0.5 with 1/v = 0.33 and (B) p = 1.0 with 1/u = 0.5. 




{/(tTV 1 /^;^). The percolation phase transitions of ran- 
dom networks under the GAP are continuous always and 
its universality class depends on the probability parame- 
ter p. 

This work is supported by the National Natural Sci- 
ence Foundation of China under grant 10835005. 



Figure 6: (7(0; p) with error bars smaller than the symbol, 
a given p and different N collapse into a scaling function 



[1] R. Albert and A.-L. Barabasi, Rev. Mod. Phys. 74, 47 
(2002). 

:] D.J Watts and S.H Strogatz, Nature 393,409 (1998) 
■] M. E. J. Newman, Networks: An introduction, Oxford 

University Press, New York (2010). 
] D. Stauffer and A. Aharony, Introduction to Percolation 

Theory (Taylor & Francis, London, 1994). 
] Bela Bollobas, Random Graphs (Cambridge University 

Press, Cambridge, 2001). 
.] S. N. Dorogovtsev, A. V. Goltsev, and J. F. F. Mendes, 

Rev. Mod. Phys. 80, 1275 (2008). 
'] D. Achlioptas, R. M. D'Sousa, and J. Spencer, Science 

323, 1453 (2009). 
,] R. M. Ziff, Phys. Rev. Lett. 103, 045701 (2009). 
'] Y. S. Cho, J. S. Kim, J. Park, B. Kahng, and D. Kim, 

Phys. Rev. Lett. 103, 135702 (2009). 
i] F. Radicchi and S. Fortunato, Phys. Rev. Lett. 

103,168701 (2009). 
] E. J. Friedman and A. S. Landsberg, Phys. Rev. Lett. 

103, 255701 (2009). 
:] Y. S. Cho, B. Kahng, and D. Kim, Phys. Rev. E. 81, 

030103 (2010). 

,] R. M. D'Souza and M. Mitzenmacher, Phys. Rev. Lett. 

104, 195702 (2010) 
: ] N. A. M. Araujo and H. J. Herrmann, Phys. Rev. Lett. 

105, 035701 (2010). 
■] S. S. Manna and A. Chatterjee, Physica A 390, 177 

(2011). 

i] Y. S. Cho, S.-W. Kim, and J. D. Noh, B. Kahng, and 

D. Kim, Phys. Rev. E 82, 042102 (2010). 
] R. M. Ziff, Phys. Rev. E. 82, 051105 (2010). 
] J. Nagler, A. Levina, and M. Timme, Nature Physics 
7,(2011) 

'] R. A. da Costa, S. N. Dorogovtsev, A. V. Goltsev, and 

J. F. F. Mendes, Phys. Rev. Lett. 105, 255701 (2010). 
i] O. Riordan and L. Warnke, Science 333, 322 (2011). 
] S. Fortunato and F. Radicchi. rarXiv:1101.3567t ^l. 
:] M.X. Liu, J.F. Fan, L.S. Li, and X.S. Chen, Eur. Phys. 

J. B, submitted. 
.] V.Privman and M.E.Fisher, Phys. Rev. B. 30 ,322(1984). 
: ] Finite Size Scaling and Numerical Simulation of Statisti- 
cal Systems, edited by V.Privman (World Scientific, Sin- 
gapore, 1990). 

■] P. Erdos and A. Renyi, Publ. Math. Inst. Hungar. Acad. 

Sci.5, 17 (1960). 
i] P. Grassberger, C. Christensen, G. Bizhani, SW. Son, 
and M. Paczuski, Phys. Rev. Lett. 106, 225701 (2011). 



[2] 
[3] 

[4] 

[5] 

[6] 

[7] 

[8] 
[9] 

[10 

[11 
[12, 
[13 

[i4; 

[15 

[16 

[17 
[18 

[19 

[20 
[21 

[22 

[23 
[24 

[25 
[26 



